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Abstract 

The development of Schramm-Loewner evolution (SLE) as the scaling limits of 
discrete models from statistical physics makes direct simulation of SLE an important 
task. The most common method, suggested by Marshall and Rohde [MR05], is to 
sample Brownian motion at discrete times, interpolate appropriately in between 
and solve explicitly the Loewner equation with this approximation. This algorithm 
always produces piecewise smooth non self-intersecting curves whereas SLE K has 
been proven to be simple for k 6 [0,4], self-touching for k 6 (4,8) and space-filling 
for k > 8. In this paper we show that this sequence of curves converges to SLE K 
for all k ^ 8 by giving a condition on deterministic driving functions to ensure the 
sup-norm convergence of simulated curves when we use this algorithm. 

1 Introduction 

The Loewner equation uses a real valued function, the Loewner driving term, to describe 
a family of decreasing simply connected domains in the complex plane. It was first in- 
troduced by Charles Loewner as an attempt to solve the Bieberbach conjecture. This 
conjecture was completely solved by de Branges with the Loewner equation as one of the 
key tools |dB85] . Oded Schramm rediscovered the Loewner equation when he was studying 
the scaling limits of discrete models. In this context, he introduced the Schramm-Loewner 
evolution (SLE K , k > 0), a random growth process in the plane |Sch00] . This process is 
obtained from the Loewner equation with a random driving term which is a/k times Brow- 
nian motion. By the work of Lawler, Schramm, Sheffield, Smirnov, Werner and others, 
SLE arises as a scaling limit of various discrete models from statistical physics [LSW04J, 
[SS05J, [SS09], |Smi01j . |Smil0j . It is therefore very desirable to generate pictures of SLE K 
directly to help understand the discrete random paths from those models. 

We are primarily interested in the case when the driving function corresponds to a 
growing curve. There are so far two methods to directly simulate the Loewner equation. 
The first method uses the fact that the Loewner equation is a first order ODE, and hence 



one can numerically solve the equation, for example using Euler's method. Some of the 
first simulations of the SLE curve were obtained by Vincent Beffara using this method. 
It produces a good approximation to the SLE K hull (not the path) for k > 4. One 
disadvantage is that it does not show the curve corresponding to the driving function but 
only a neighborhood of it, see Section [2^2 We note that SLE K is a random curve for all 



k > 0, see [RS05] and [LSW04J. This approach has not been often used, and we do not 
discuss its convergence here. 

The second method for simulating SLE was suggested by Marshall and Rohde [MR05J . 
The algorithm discretizes the driving function and square-root-interpolates it. As a result 
the algorithm approximates SLE maps by composing many basic conformal maps, which 
are easy to compute. The algorithm was described and implemented in |Ken07j . |Ken09] 
as well as in many other works. Oded Schramm was skeptical at first that the pictures 
generated from this algorithm well-present the SLE curves, according to Steffen Rohde 
|Rohj . The curves simulated from the algorithm are piecewise smooth and simple, see 
Figure [TJ whereas SLE K is a random fractal curve which is simple for k G [0,4], self- 
touching for k G (4, 8) and plane-filling for k > 8. In this paper, we will prove: 

Theorem 1.1. For k / 8, let 7™ be the sequence of curves simulated from the second 
algorithm. Then under the half plane parametrization, the sequence 7™ almost surely 
converges to SLE K in the sup-norm. 

It is known that, for all k, the sequence 7 n converges to SLE in the context of 
Caratheodory convergence |Law05j and Cauchy transforms of probability measures |Bau03j . 
However these types of convergence relate to Loewner chains rather than curves, see 
|Law05| Chapter 4] and respectively |Bau03j for details. For k < 4, when one views 
curves as compact sets, the sequence 7" converges almost surely to SLE K in Hausdorff 
metric |BJK12} Section 7}. A general principle is to set up a theorem for the determin- 
istic Loewner equation and then translate the result into the SLE context. The theorem 



1.1 will follow from a more general theorem for deterministic curves. In particular, we 



will show that there is a class of driving functions for which the sup-norm convergence 



of approximation curves occurs, see Theorem |2.2| for the details of the statement. It is 
shown in |MR05| . |Lin05j and |LR10] that driving functions whose H61der-l/2 norms are 
less than 4 generate simple curves and that the Hilbert space filling curve is generated by 



a H61der-l/2 function. Our Theorem 2.2 is also applied to these driving functions. 



Corollary 1.2. Consider the driving function that generates the Hilbert space filling curve 
or that has Holder- 1 /2 norm less than 4- Then the sequence of curves simulated from the 
algorithm for this driving function converges uniformly. 



We note that Theorem 2.2 also provides the convergence rate of the simulation. The 
key is to estimate how the curve changes when we modify the driving function since the 
driving functions of simulated curves converge uniformly. There are two key estimates 



in the proof of Theorem One is the boundary behavior of a conformal map and the 
other is the perturbation of a Loewner chain when there is a small change of its driving 
function. The latter is a Gronwall-type estimate which appears in |JV12j and |JRW12] . 
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The two estimates are both related to the growth of the derivative of conformal maps 



near the boundary which will provide to the assumptions of the main theorem 2.2 



The paper is organized as follows. In section [2j we begin by reviewing the Loewner 
equation; we then state our main result and prove some preliminary lemmas. In section 



2.2, we describe the first and second algorithms simulating the Loewner equation. Then 



the main theorem is proved in section |3j The applications will be discussed in Section [4] 
when we show the convergence rate and discuss several variants of the algorithm. 

Acknowledgment. The author would like to thank Steffen Rohde for numerous in- 
sightful conversations. The author also would like to thank Elliot Paquette and Brent 
Werness for helpful comments on early drafts of this paper. 
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Figure 1: Simulations of SLE 8 /3 (left) and SLE 6 (right) from the same Brownian motion 
sample with 12800 points 



2 Loewner equation, algorithms, main result 
2.1 Loewner equation 

There are many versions of the Loewner equation. In this paper we focus on the (down- 
ward) chordal Loewner equation in the upper half plane H: 



d t gt{z) 



9t(z) - A(t) 
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with the initial condition go(z) = z for every zGH, where A(i) is a real- valued continuous 
function defined for t > 0. Sometimes we write A^ for A(t). The family of (gt)t>o is called 
the Loewner chain. 

For each point z G H, the solution of ([I]): t 4 ^(z) is uniquely defined up to 
T 2 = inf{t > : Qt{z) = A(t)}. As t increases, the set K t = {z G EI : T z < t}, called 
the hull, grows. It is known that for each t > 0, gt is the unique conformal map from 
H t := M\K t — > EI satisfying the hydrodynamic normalization at oo, 

lim [g t (z) - z\ = 0. 
It is usually easier to work with the upward Loewner equation: 

d t h t {z) = H ~ 2 _ m , h(z) = z, (2) 

for z G EI and real- valued continuous function 

It is not hard to show that if g t , < t < T is the solution to Q with a driving 
function A and h t is the solution to ^ with = A(T — t) then 

h T {z) = g T l {z). 

We are interested in the case that the Loewner chain (g t ) is generated by a curve 7, 
i.e., H t is the unbounded component of HP\7([0, t\). It follows from Theorem 4.1 (JRS05J) 
that this is equivalent to the existence and the continuity in t > of 

P(t) := lim . g;\\{t)+iy). 
2.2 Algorithms simulating Loewner equations 

Let us briefly discuss the first algorithm mentioned in the introduction. This idea to 
simulate K t is to determine whether a point z in the upper half plane EI satisfies T z < 
t. One cannot examine all the points in EI so if T z < t one declares that a certain 
neighborhood of z is in K t . To calculate the blow-up time T z one needs to run equation 
until gt(z) hits A(t). However, there is no general method to solve with given 
driving function. There are a few cases one can solve explicitly, see |KNK04j . As a result, 
if 7 is the simple curve corresponding to A then the simulation of 7([0,t]) = Kt is a 
neighborhood of the actual 7QO, t}). This is often not a good way to visualize the curve 
7- 

We now discuss the second algorithm to simulate Loewner curves. It was first ap- 
peared in |MR05j . The algorithm has also been described in |Ken07j . |Ken09j . where 
modifications and fast implementations are discussed. One advantage of this algorithm is 
that it always produces simple curves. For the rest of the paper, this is the algorithm we 
consider, unless otherwise stated. 
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The algorithm is based on two observations. First, fix s > 0, and let (gt) be the 
solution of the Loewner equation with driving function X(t) = X(s +t),t > 0. This 
solution can be obtained by g s+t o g^ 1 . Indeed 

dt9s+t° g^(z) 



9s+t o g l {z) - X(s + t) 9s+t o g-i( z ) - X(t) ' 

and g s ogj x (z) = z - By the uniqueness of solution of the equation ([T|, gt(z) = g s+t ogg 1 (z). 
If we let K t be the hull associated with g t then 

9s(K s+t ) = K t and K s+t = K s U gj\k t ). 

So in order to compute K s+t , one can compute K s and g" 1 , by using the information of 
A on [0, s], and compute K t by using A on [s, s + t\. 

The second observation is that when A is of the form c\ft + d, for some real constants 
c and d, one can solve for K t explicitly. In this case, K t is a segment in the upper half 
plane starting at d G R that makes an angle air with the positive real axis where 

1 1 c 

a — , 

2 2 V /16+^' 

and gt\z + A(t)) = (z + 2 v / i y ^^) 1 - a ( y z - 2^t^~^) a + d. See |KJNK04j for a proof. 



We now fix a step n > 1. Let = | for < k < n. So to — 0, ti, ■ ■ ■ , t n — 1 is a 
partition of [0, 1]. We will solve the Loewner equation with driving functions X(t + t k ) 
for < t < -. By the remarks above, one should approximate these driving functions 
by cy/i + d so that one can solve explicitly the Loewner equation. More specifically, we 
approximate A by A n such that they attain the same values at tk's and that A n is a scaling 
and translation of \/t on ifc+i] from A(t^) to A(tfc +1 ). Hence the function A n is defined 
as follows: 



X n (t) = y/n(X(t k+1 ) - X(t k ))Vt -t k + A(tfc) on [t k ,t k+1 ]. 
This driving function always produces a simple curve 7" : [0,1] ->HU {A(0)}. 

Denote by (<7")o<t<i the Loewner chain corresponding to A™. Let /" be the inverse 
function of g[ l and ft(z) = f™(z + X n (t)). Define 

so that 

/tfc = G fc-i oG fc-2 o---oG . 
For each t G [0, 1], let 7" be the image of 7" under g™ — A n (t), i.e., 

7r(^)=^(7 n (t + ^))-A n (t) and 7 »(t + a) = ^ (7^)) for < 5 < 1 - t. (3) 

We have chosen A n so that 7™ ([0, -]) is a segment starting at and that G k has an explicit 
formula: 

l-a 



G n k (z)= I ; • 2, ' " ' ' ; " 



na ) \ y n(l — a) 



n 

9t k 



A n (0) 




\ n (t k ) 



z + X n (t k ) 



'tfc + i S 











Figure 2: At each step k, we compute G%, / t " and 7?. The fc-th sub-arc of the simulation 
curve 7 n is the image of 7^ under / t ™ . 



where a 



1 _ 1 yH(A(f fc+ i)-A(t fc )) 

2 2 A /l6+n(A(i fc+1 )-A(t fc )) 2 



e (0, 1). See Figure 



Therefore in order to compute 7 n ([0, 1]), we find 7?; ([0, -]), f™, and then 

7"W = ft(lt(t - t k )) for t E [t k ,t k+1 ),0< k < n - 1. (4) 
Notice that A n converges uniformly to A on [0, 1] since 

sup |A n (t) - \{t)\ < 2 sup |A(t)-A(s)|. 

te[o,i] s ,ts[o,i],|t-s|<i 

We mention without proof a geometric property of GjJ which we will use later. 

Lemma 2.1. Consider the conformal map G(z) — {z + a) 1-Q (z — 6) a /rom EI to EI minus 
a slit starting at 0, where a, b > 0, a G (0, 1) and aa — (1 — TTie point is mapped to 
the tip of the slit. Then the imaginary part of G(iy) is increasing on (0, 00). In particular, 
the image of iy has a larger imaginary part than that of the tip of the slit. 



The dashed line of Figure [2] illustrates this proposition. 
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2.3 Main results 



We shall consider driving functions which have the same regularity as Brownian motion. 
A subpower function is a non decreasing function from [0, oo) to [0, oo) satisfying: 

lim x~ u <f)(x) = for all v > 0. 

If </>i, </>2 are subpower functions then so are ctfix, <p\ and max(0i, 02 ) for every c > 0. 

The function A is called weakly Holder- 1/2 if there exists a subpower function tp such 
that 

osc(A;5) := sup{|A(t) - \{s)\ : s,t G [0,1], |* - s| < 8} < V5(p(l/8) for all 8 > 0. (5) 

It follows from P.Levy's theorem that the sample paths of Brownian motion are almost 
surely weakly H61der-l/2 with subpower function c-y/log(<5), c > y/2, see |RY99t Theorem 
1.2.7]. 

It is known that if A is weakly H61der-l/2 and if there exist cq > 0, yo > and 
< /3 < 1 so that 

\f' t {iy)\ < c y~P for all < y < y , t E [0, 1], (6) 



where f t {-) = ^ _1 (A(t) + •), then (gt)o<t<i is generated by a curve, see |JVLllt section 
3]. This is one of the main ideas to show the existence of SLE curves for QRS05J). 
We note that the Loewner chain of SLE 8 does not satisfy (J6]). 

Our main theorem shows that under these hypotheses the algorithm gives the sup- 
norm convergence of the simulation curves. 

Theorem 2.2. Suppose X is a weakly Hdlder-1/2 driving function with a subpower func- 
tion cp and suppose the condition is satisfied. Then the curve 7 generated from the 
Loewner equation can be approximated by the algorithm; that is, there exists a subpower 
function <p such that for all n > ^ and t e [0, 1], 

l7»W-7«l<-f^, (7) 

""■2 I 1 V 2 I 



H2' 



where 7™ is the curve generated from the algorithm which is explained in section 2.2 The 
function cp depends on ip, c and (3. 



A related question is the following: under what additional assumptions, does the 
uniform convergence of driving functions imply convergence of corresponding curves? A 
priori the convergence of curves occurs in the sense of Caratheodory convergence, see 
|Law05j . and in the sense of Cauchy transform of probability measures, see [Bau03j for 
definitions and details. As said in the introduction, these types of convergence do not 
directly involve to the curves. In |BJK12l Section 7], it is shown that if two driving 
functions generating simple curves are close in the sup-norm and if one function has the 
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condition ^ then the two generated curves are close in Hausdorff distance. One really 
wants to see two curves are close in the sup-norm. Lind, Marshall and Rohde [LMR10] 
show that if the driving functions have H61der-l/2 norm less than 4, then the curves 
converge uniformly. However, the Brownian motion is a.s. not Holder-1/2. In |SS12j . 
the authors study sufficient conditions to have uniform convergence of bidirectional paths 
(the curves and their time-reversals). The paper by Johansson- Viklund |JV12] uses the 
tip structure modulus to get another criterion for uniform convergence of curves. 

In the rest of this paper, C stands for absolute constant and for general subpower 
functions; c and ip stand for constants and subpower functions that may depend on the 



assumptions of Theorem 2.2 They can change line by line and are indexed when necessary 
to avoid confusion. 

Since we are interested in the same type of driving functions as those in |JVLlll 
section 3], there are several results from their paper we will use and state here for the 
convenience of the reader. 

Lemma 2.3. ]JVL11[ Lemma 3.4] Let K be a hull. There exists a constant C < oo such 
that 

hcap(K) < C diam(K) height(K) . 

Lemma 2.4. \Wonl2\ Lemma 3.1] Suppose 7 is the curve generated by a driving function 
X(t) in (7p. Then for all z E j{[0, t)), 

\Rez\< sup |A(r)-A(s)| 

0<s<r<t 

and 

Im z< 2\ft. 

Lemma 2.5. \JVL11\ Proposition 3.8] Let (g t ) be the Loewner chain corresponding to 
A(t) satisfying and |#|). Then there exists a subpower function (fx such that ifO<t< 
t + s < 1 and s e [0, y 2 } 



where 



and that 



\-y(t + s) - 7 (t)| < (p 1 (l/y)[v(t + s,y)+v(t,y)], 



v(t,y) = J* \fi(ir)\dr < ^y 1 '^ °<y<yo, 

| 7 (t + s)- 7 (f) I < ip^l/y)^-^ 1 -? for0<s<y 2 < y%. 



Most of the time we will deal with the behavior of conformal maps near the real and 
imaginary lines. For every subpower function <p, constant c > and n6l, define 

4>(n) ,1 c 
Ai,c,4 = {x + %y e M : \x < and < y < 
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Lemma 2.6. There exist constants a > and d > such that if Z\ and z 2 are inside the 
box An^, and f is a conformal map on HI then 

\f{z x )\<d<j>{nr\f{ilmz x )\ (9) 

and 

dw,hy P {zi, z 2 ) < d log0(n) + d. 

The constants a and d depend only on c, not on <p or n. The notation dn,hyp{zi, z 2 ) means 
the hyperbolic distance between z\ and z 2 in a simply connected domain H . 

Proof. The proof is similar to |JVL114 Lemma 3.2]. □ 

Lemma 2.7. \Pom9^, Corollary 1.5] (half plane version) If f is a conformal map ofM 
into C and if zo,z\ G EI then 

l/Oi) - f( z 2)\ < 2\(Imz 1 )f(z 1 )\exp(Adu,hyp(z 1 ,z 2 )). 



3 Proof of main theorem 2.2 



3.1 Heuristic argument 

Let 74 be the image of 7 under g t — \(t), i.e. 

lt(s) = g t (l(t + s)) - X(t) and j(t + s) = f t (lt(s)). 
We want to estimate 

I7OW1) - l n (h +1 )\ = \f tk (z) - fl(w)\ < \f tk (z) - f tk (w)\ + |/ t » - f»(w)\ (10) 

where z = 7 tfc (i) and w = 7£(J)- 

First, z and w are the tips of two curves generated respectively by two driving functions 
defined on [0, -]. It follows from Lemma 



2.4 



that Imz and Imw < -7=. 

— V" 



For the first term in the RHS of (10), it follows from Lemma |27T that 
\k( z ) ~ k( w )\ < 2\(Imz)fl k (z)\exp(4d u>hyp (z,w)). 

Notice that z and w trivially have positive imaginary parts. 

If we can show that z and w are in the same box A nfi ^ then combining with a 
hypothesis of f[ on zIR + , we obtain similar inequalities for \fL(z)\ and \f' t (w)\ from 
Lemma [2.61 Then it follows that 



2lmz \fl k (z)\exp(4:d Mthyp (z,w)) < (Im^) 1 "^^^) < - _ p -> as n -> 00, 
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where the notation / < g means that / < Cg for some constant C > 0. 

Since w is a tip of a straight line generated by a nice driving function, we can show 
to is in a box A nfi ^ see (13). However, z = 7t fc (-), in the case of SLE, has continuous 



density on the strip {x + iy : x G M, < y < 2}. So there might not exist a controllable- 
sized box v4 n c that contains z. However Lemma |3.2| shows the existence of a point in 



7t fe ([0, -]) PI A n ^, and we will use this point instead of Jt k {-)- Then we use the uniform 



continuity of 7 to get back to (10). 



For the second term in the RHS of (10), notice that 



/*» - /r» = /*> + a(**)) - /r> + A(t fc )). 

This expression is a perturbation of two solutions of the upward Loewner equation ([T]) 
with two driving functions 1 1— > A(tfc — •) and t y A n (t& — •). 



We will use the following lemma from [JRW12] and the fact that A and A n are close 
on [0, £*;] and that \f' t (w + \(tk))\ is well-controlled. 



Lemma 3.1. URW12L Lemma 2.3] Let < T < 00. Suppose that for t e [0,T], / t (1) and 
ff® satisfy the upward Loewner equation ^2) with and respectively, as driving 

terms. Suppose that 

e= sup \WP - WP\. 



se[o,T] 



Then if u = x + iy e HI, then 



If 



u 



f?\u)\ < eexp 



b fr,l(ffl)'(«)l ^IW, 



l;'-f( 2 ) 



1/2 



+ log log 



l T,y 



where Ir, y = a/4T + y 2 . 

Thus if K/^^HI < cy- /3 , then 

If one can show furthermore 

4>(n 



e < 



and y = Imu = Imw > — — 
n 6(n)\/n 



then 



\f?\u)-f¥(u)\ 



(2), 



0(n) c 



< 

V 2 > 



0. 



From here, we only have an estimate for Yi{t) — 7 n (£)| when t = tk- To have an estimate 
on the whole interval we notice that 

7"([w fc+2 ]) = fZKiK 2 -\) = gk^jo, 



n n 
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It follows from a property of G% (Lemma 2.1), that every point in 7t!([-, f]) is in 
a box A nfi ^ and hence we can apply the same argument for (10) with 7 n (tfc + i) being 
replaced by any j n (t), t k+1 <t< t k+2 . 

Now we will go into the details of the proof. 



3.2 Proof of Theorem [2^2 



Fix an arbitrary interval / = [t k , tk+2], < k < n — 2. Denote 7^ = 7t fc , 7^ = 7^. 

We will estimate \^{s + t k ) — 7 n (r + t k )\ for all r G [-, -] and with a specific s chosen 
later. Combining with the uniform continuity of 7, we will have an estimate for 

1 2, 



h(r + t k ) - 7 n (r + £ fc ) I with all r G [ 



n n 



From now on, we will choose n so that - < y$. Denote z = 7fc(s), w = 7^(r). By the 
triangle inequality, 

\l(s + h) - 7 "(r + t k )\ = \f tk {z) - t(w)\ < \f\ k (z) - / t »| + |/ t » - t{w)\. (11) 



The first term in the RHS of (11). It follows from Lemma 2.7 that 



\f tk (z) - f tk (w)\ < (2lmz)\n k (z)\exp(4d u , hyp (z } w)). (12) 

The next lemma shows the existence of a point in 7fe([0, -]) D A nfi ^. 
Lemma 3.2. There exists a subpower function depending only on ip, Cq and (3 of The- 



orem 



~2J^ such that for n > 1 and < k < n — 1, there exists s G [0, |] such that 



lk{s) G 

Proof. Since rj := 7fe([0, |]) is the curve generated by the Loewner equation (jlj) with 



driving function X(t k + .) — X(t k ) on [0, -] and since A is weakly H61der-l/2, it follows 
from Lemma [2.41 that 



Re7 fc s < \ -cp - =: —=r- and Im 7A; (s < —= 
V n 2 Jn \/n 



for all s G [0, f]. 
This implies 



, . , 1 / \ , w n 2^ 2(p 2 (n) ipJn) 
diam{rj) < height (77) + width(?7) < — + 



jn \/n \/n 

It follows from Lemma [2.31 that 

2 

— = hcap(?]) < C diam(r^) height(?7), 

n 
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so, 



height (77) > 



The lemma follows by choosing a highest point and letting <j) := ^5 := max(<^4, ipz). □ 



With this specific point 7&(s), we can use the inequality ^ in Lemma 2.6 To have a 
bound for exp(4dm,hy P (z, w)) one needs to show that w = 7^ (r) is also in A n% ^ . By the 
same argument in the above lemma, since J k ([0, ^]) and 7fc +1 ([0, ^]) are line segments, 



and 



g A 



(13) 



Lemma 3.3. There exists a subpower function <fi depending only on ip, Cq and (3 such that 
for all n, k and r G 7^(r) zs in the box A„ >2 ^. 

Proof. Notice that 7^ (r) is the tip of the Loewner curve generated by the driving function 

(14) 



t ^ \ n (t + t k ) - X n (t k ), t G [0, r). It follows from Lemma 2.4 that 



|Re 7fc »| <suv{\\ n (t + t k )-\ n (t k )\,te [0,r]} < J^(^ 



n 



and 



Im7fc(r) < 2Vr < 2 



2.6 



Now the rest of the proof is to find a lower bound for Im^(r). Fix r G [-,-]• Let 
x + iy := ! n k+l {r ~ £), G := (/fj" 1 o 

Since 7^ +1 ([0, ^]) is a line segment with the tip 7^ +1 (l/n) in A n2 ^2 , by Lemma 

d M>hyp (x + iy,iy) < C dog <p 5 (ri) + C 

and 

d m ,hy P (x + iy,iy) = G? HVr « [0 1^(7^), G(m/)) > d Mjhyp (^(r),G(iy)). 



It follows from (13) and Lemma 2.1 that 



Im 7fe » > lmG{W l > lmrk{lM > 



C(f 5 (n) c C(f 5 (n) c Cy/n(p 5 (n) c+1 y/ntp 6 (n)' 
So 7^(r) and j k (s) are both in A n ^ m with (p 7 = max(y2 5 , <p 6 ). 



□ 



We now apply Lemma 2J3 and obtain 

\ftk(z) ~ ft„H\ < {21mz)\fl k {z)\exp{4d mMp {z,w)) 

< C(Imz)ip 7 (nT\fl k (rtmz)\exp(C\og<p 7 (n) + C) 

< (Im z) 1 ^ Cc Q if 7 (n) a exp(C log if 7 (n) + C) 



< 



n 



-PCco<p 7 {n) a exp{C\og<p 7 (n) + C) 



1-/3 ' 



(15) 
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The second term in the RHS of (11) 



Let u = x + iy := w + A(tjt). Since A(ifc) = \ n (tk), 

4H - fl(w) = f tk (u) - /»(«). 



Applying Lemma 3.1, we get 



!/**(«) -/£(«)! <e«p 



1 



log 



log 



1/2 



+ log log 



where I tkfV = ^ At k + y 2 and £ = sup te[0)tfc] \\{t) - X n (t)\ < ^gi 
Since y = Imw = Imw G [ — - — 



y/rup-tin)' 1 -^/jiJ' 



y 



mp 7 {n). 



Since / t ' fc (w) = /^H, 



Also 



Co 



K/rj'WI < ^(r 1 + 1) < 2C<p 7 {n)Vii, 



where the first inequality comes holds for all hydrodynamic normalized conformal map of 
EL 

It follows that 



\kH-t(w)\ < ^exp 



n 



1 + 13 



\og(c(f7(n)y/n) + log log 2 v Imp^in 



n 



2 



(16) 



End of the proof of Theorem 2.2 



It follows from (11), (15) and p§ that 



\l{s + t k ) - 7 n (r + t fc )| < 



■1-/8 



+ 



<p 9 (n) 



<Pio(n) 



n) 



l- 



n) 



1-1 



for all re [± 

Using the uniform continuity ^ of 7, we obtain 



|7(r)-7"(r)|< 



1 , n+w 
n) V 2 



for all r G [tfc+i , tfc+2] an d < < n — 2, hence for all r G [0, 1]. 



□ 



Proof of Corollary |1.2 



It was shown in [MR05J, [Lin05] and |LR10] that the unbounded complements of the hulls 
generated by the driving function are John domains. Therefore, it follows from |Pom92[ 
Chapter 5] that the condition (fsj) is satisfied. □ 
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4 Applications 



4.1 Variants of the algorithm 



Variant 1. The conclusion of Theorem |2.2| still holds for every A n that satisfies: 

(p{n) 



\X n (t k ) - X(t k )\ < 



n 



(17) 



and 



X n (t) = yfn{\ n {t k ) - X n (t k ))Vt^T k + \ n (t k ) on [t k ,t k+1 ]. (18) 

Indeed, the main inequality pl| has a slightly change 

| 7 (a + t k ) - 7»(r + t*)| = |/ tfc (*) - < \f tk (z) - f tk (w + X n (t k ) - X(t k ))\ 

+ \f tk (w + X n (t k ))-f? k (w + X n (t k ))\. (19) 

We can see that z and w + X n {t k ) — X{t k ) are still in the same box A nfi ^. Hence the 
same argument follows. 



Variant 2. Lemma 2.1 and the property (13), hence Lemma 3.3 are still true if 



instead of square-root-interpolating X n on [t k , t k +i], we consider any function interpolating 
between X n (t k ) and X n (t k+ i) such that when we run the Loewner equation for < t < -, 
the resulting 7^ has non decreasing imaginary part on [0, -]. 

In particular, linear-interpolating X n 

X n (t) = X n (t k ) + n(X n (t k+1 ) - X n (t k ))(t - t k ) on [t k , t k+1 ] 



also give the same conclusion as in Theorem 2.2 (see |KNK04] for linear driving functions). 



Variant 3. Instead of using tilted slits on each small interval one can use vertical 
slits [Ken09j. In this case, A™ is a step function: 

X n (t) = X(t k ) for t G Mk+i) 



and 



It 



?(t) = X(t k ) + 2iVion [0,-). 

k n 



However 7" defined by Q is not a curve. We can do as follows. We compute 7™ at discrete 
points t — to, ti, ■ ■ ■ , t n . Then connect them with a straight line in that order. As plotted 
in |Ken09] . it is almost impossible to distinguish this curve and the one from the (main) 



algorithm. Indeed, the same proof of Theorem 2^ is carried over for this algorithm and 
the same conclusion of this theorem holds. 
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4.2 Speed of convergence to SLE n 



We can estimate the speed of convergence of the algorithm to 7 K := SLE K with k 7^ 
Corollary 4.1. There exist constants C\, 02,03,04 > depending on k such that 



Ci(logm* 



f'2 



P I |7 K - 7 m | l[o,i],oo < v °_ ' /or a« m > n > 1 
\ x/m 1 ^ 2 / 



c 3 



71 



In other words, this corollary implies Theorem 



We will apply Theorem 2.2 to the case X(t) = \fnB t . It follows from |LL10| Theorem 
3.2.4] that there exist constants c\ (depending on k) and c 2 such that 



m 



P l osc(A; — ) > Ci\/ — for all m > n \ < 



m 



rr 



(20) 



Notice that in Theorem 2.2, if ip(n) = y/\ogn in (nj) then goes through the proof, the 



subpower function in (M) is of the form c(logn) c ' for some constants c and d. 



It follows from |JVL111 Proposition 4.2] that there exist constants 0' G (0, 1), C3 and 
C4 > depending on k such that 



00 2 Z 



2^2^ p {l/o--i)2--(^ m )l> 2 " m )< 

m=n j=l 



_C3_ 

2«C4 



or 



This implies that 

P{|/y_i )2 - 2 m(i2- m )| < 2 /3 ' m for all 1 < j < 2 2rn ,m> n} > 1 

Thus there exist C5 > and (3 G (/?', 1) such that 

P{|/;(^)| < ^ for all < y < 2~ n , t G [0,1]} > 1 - 



P<| \fl(iy)\ < % for all < y < -±=,t G [0,1] |> > 1 ' ' 



_C3_ 



y 



n 



TV 



Ci/2- 



Combining (20), (21) and Theorem 2.2, we get 



P lb" - 7™ ll M « < ^£ for a., m > „ > 1 - (| + ^) 



(21) 



which proves Corollary 4.1 



□ 
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4.3 Random walk algorithm to simulate SLE curves 

This algorithm |Ken09l Section 2] is based on the Donsker's invariance theorem: a scaling 
limit of simple random walk converges in distribution to the Brownian motion. 

For fix k > 0. We choose a 6 (0, |] such that 

4(1 -2a) 2 

K = — . 

a(l — a) 

Let fi(z) = (z + l — a) 1 ~ a (z — a) a , f2(z) = (z + a) a (z — (1 — a)) 1- ®. For every i > 1, choose 
4>i = fi or (pi = / 2 with equal probability. Then we compute inductively F n = F n _i o n 
with Fq = id. The map F n is conformal from EI to EI minus a slit curve. After rescaling 
and translating so that this slit curve has the half plane capacity 1, we get a simple curve 
7™. More explicitly, 7 n is generated by A n whose formula is 

S 

X n (t k ) = s/k^ for all t k , 
in 



and \ n {t) = y/n(X n (t k+ i) - \ n (t k ))Vt zr t k ~ + A n (t fc ) on [t k ,t k+1 ], 
where S k = X x + • • • + X k , X[s are iid and PpQ = 1) = P(X< = -1) = \. 

By Donsker's invariance theorem, A n -4 i/kBI^ on C([0, 1], ||.||oo) • So Et\7 n ([0, 1]) -4- 
H\7 K ([0, 1]) in the context of Caratheodory kernel convergence |Law05] and Cauchy trans- 
forms of probability measures [Bau03j. Kennedy |Ken09j raised a question whether 7 n 
converges in distribution to r y K . 

We now show that 7 n converges in distribution to 7 K under the sup-norm of C([0, 1]) 
when k 7^ 8. Indeed, it follows from |LL101 Theorem 7.1.1] that for each n, we can couple 
A" and the Brownian motion in the same probability space such that 

P{ max | X n (tj ) - yfcB t . \ > C ^^ gH } < Cn~ 3 (22) 
for some universal constant C > 0. 



Hence from (20), (21), (22) and the discussion of Variant 1, there exist constants c% 
and eg depending on k such that 

-m ^ c 8 (logn) C9 [ ^ c 2 c 3 C 
MIL' -7 ll[o,i],oo< ^/^p 1 -^-^-^- 



This implies that 7™ converges in distribution to 7 K . 
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